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Abstract 



We recall the main features of the Regge approach used to understand soft interactions 
at LHC and higher energies. Unitarity tames the power growth of the elastic proton- 
proton scattering amplitude with energy, and leads to the migration of the secondary 
particles produced in high-energy proton-proton collisions to larger transverse momenta. 
We discuss, in qualitative terms, the role of processes containing large rapidity gaps 
^ ■ (LRG), and the probability that the gaps survive population by secondaries produced in 

additional soft interactions. We explain how the Regge diagram corresponding to a LRG 
event simultaneously describes events with different (single, double, etc.) particle density 
in the same rapidity interval. We show that the role of these, enhanced, multi-Pomeron 
diagrams can be studied by measuring multiplicity fluctuations and long-range rapidity 
correlations between secondaries produced at the Tevatron and the LHC. Finally, we make 
a list of the characteristic features of the multi-Pomeron description of soft interactions 
that may be observed at the high energies accessible at the Tevatron and the LHC. 



1 Topical Review to be published in J.Phys.G, based on two lectures by Misha Ryskin at the St. Petersburg 
Winter School, Gatchina, February 2009 



1 Introduction 



"Soft" physics is manifest in various places in different guises. The interest in soft interactions 
at very high energies is stimulated, theoretically, by the fact that the observed hadronic cross 
sections are growing with energy and, experimentally, by the advent of the LHC. 

At the LHC energy we have to allow for unitarity effects, which are necessary to provide 
the consistency of the strong interaction. In fact, when we face a cross section which grows as 
a power of the energy, there are two possible ways to restore unitarity. 

First, we may introduce a new particle (sufficiently heavy to have not been observed yet) 
which cancels the growth of the cross section at energies larger than the mass of the particle 
(for example, Higgs boson, super-partners, etc.). After this, we deal with a weak interaction, 
where the main contribution comes from one, or a few, simple diagrams. More complicated 
Feynman graphs just describe corrections to the lowest-order amplitude. 

Another possibility is to consider a strong interaction. Here, from the beginning, we start 
with a Hermitian Lagrangian, which should already account for unitarity. When the Born 
amplitude (or an amplitude obtained via the summation of some group of diagrams) becomes 
too large and the interaction becomes strong, then new more complicated diagrams must enter 
the game. It is the contribution due to these new graphs which tames the growth of cross 
section such that the final result satisfies unitarity. 

In the case of QCD, where the BFKL amplitude (A oc s 1+A ) grows as a non-integer power 
(A oc a s ) of energy, y/s, we expect to observe the second scenarico Therefore it is important to 
observe, and to study, the role of more complicated diagrams (arising from multiple interactions) 
at LHC energies, and to trace how the theory restores the unitarity of the strong interaction. 
That is, from a theoretical viewpoint, it is of great interest to observe experimentally how 
unitarity tames the growth of high-energy hadronic amplitudes, leading to saturation both 
in the transverse momenta k t and in the impact parameter b distributions of the produced 
particles. In other words, unitarity replaces the growth of the amplitude by the growth of 
configuration space occupied by the particles, both in b and in k t . 

Here, we do not present new results but attempt to recall, using modern language, the ideas 
and understanding of high-energy soft interactions which originated some 40 or more years ago. 
The discussion may be divided into three main topics: (i) elastic scattering and the total cross 
section, (ii) the cross sections of processes with gaps in rapidity (including high-mass diffractive 
dissociation), and, (iii) the saturation of particle densities in transverse momenta kt- 

Soft high-energy pp interactions are clearly important at the LHC. Moreover, they can 
complicate our ability to observe new physics. For example, topic (iii) is relevant to the search 
for new physics, since particles with a rather large k t from the underlying inclusive event will 
affect the jet searching algorithm, as is already the case for jets at HERA and the Tevatron. 
On the other hand, the evaluation of the cross section for an exclusive process, such as pp — > 

2 It looks impossible to cancel the non-integer power of s by introducing a new particle. 
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p + Higgs + p, requires knowledge of the gap survival probability S 2 , which is the subject of 
topic (ii), see, for example, [TJ [21 El H]- 

Here, we focus on the qualitative features of high energy soft interactions. These features 
are quite general. With the advent of the LHC, it is timely to gather them together and to 
emphasize the underlying physics. 



2 Total and elastic cross sections 

The behaviour of scattering amplitudes in the high energy, y/s, small momentum-transfer 
squared, —t, domain is well described by Regge theory; that is, by the singularities of am- 
plitudes in the complex angular momentum, j, plane, see, for example, [5JE]. For instance, the 
measured ir~p — > TT°n amplitude behaves as 

A{s,t) oc s a '® (1) 

where the p-trajectory, j = a p (t) ~ 0.5 + 0.9t (with t in GeV 2 ), passes through the spin-1 
p-meson resonance in the 'crossed' t-channel 7r~7r° — > pn; that is, a p (t = m 2 p ) = 1. 

On the other hand, total cross sections are observed to grow slowly with energy and are 
associated with the exchange of a trajectory with vacuum quantum numbers. The simplest 
possibility is to assume that at high energy these cross sections, such as the pp total cross 
section er tot , are driven by an isolated pole at j = a(t), which gives an (pp) elastic amplitude 

A(s,t) oc s ap{t \ (2) 

and, via the optical theorem of Fig. CD, a total cross section 

a tot oc s a ^~\ (3) 

The pole with the largest intercept, originally with ap(0) = 1 since high energy total cross 
sections were thought to have constant asymptotic behaviour, was called the PomerorlE Here 
we are interested in Tevatron and higher energies which are sufficiently large to be able to 
neglect the contributions of all secondary trajectories (which all have intercepts a(0) ~ 0.5). 

A popular parametrization of elastic pp-scattering amplitude by Donnachie-Landshoff (DL) 
is the Regge form [7] 

A(s,t)= Wo F 2 (t)s a ^ (4) 

where the signature factor r] gives the complex phase, F\ is the electromagnetic form factor of 
the proton, and the effective soft Pomeron trajectory 

a P = 1 + A + a't ~ 1 + 0.08 + 0.25t, (5) 



3 A discussion of the history of the Pomeron is given in [6]. 
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Figure 1: (a) A sketch of the optical theorem which, via unitarity, relates the pp total cross 
section to the imaginary part of the forward elastic pp amplitude; (b) the application of the 
optical theorem to obtain the cross section for the dissociation of the 'beam' proton into a 
high-mass system, see Sections |3] and 15. 11 The high-energy Pomeron-exchange expressions that 
are shown for the bare amplitudes have sizeable absorptive (that is, unitarity or screening) 
corrections; g^ and g^p are the Pomeron-proton and triple- Pomeron couplings respectively. 



where t is in GeV 2 . The intercept cn(0) just above 1 reproduces the observed slow growth 
of the total hadron-hadron cross sections at high energies. Indeed, the effective Pomeron pole 
amplitude, (jl]), gives a good description of the total and elastic differential cross section data up 
to Tevatron energies. However, this simple parametrization is expected to become increasingly 
deficient at higher energies. Unitarity is an easy way to see this. 

To discuss unitarity effects it is convenient to work in impact parameter, b, space, since at 
high energy the position of the particle in b is practically frozen. Thus the value of b determines 
the orbital momentum I of the incoming proton, I = y/sb/2; that is, fixed b corresponds to a 
particular partial wave I. In b space the elastic unitarity equation^ 

2ImA(6) = \A(b)\ 2 + GW(&) (6) 

limits the value of Iul4 < 2. However, this limit corresponds to a pure elastic interaction with 
phase of elastic amplitude, 

^ = z(l-exp(2^)), (7) 

equal to 5i = tt/2. Normally at high energies the inelastic contribution Gi nc i dominates, leading 
to a large inelasticity; that is, to a large imaginary part, Im 5 > 1. In this so-called 'black 

4 The amplitude is normalised to <7tot = 2 J d 2 b lmA(b). 
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disk' limito, we have ImA — > 1. Recall that, if we were to neglect the small contribution from 
the real part, then 1mA (b) can be measured directly from experiment as the Fourier transform 

where g 2 = \t\, Jq is a Bessel function and p 2 = (KeA/lmA) 2 < 0.02. 

Note that the DL amplitude A(b = 0) crosses the black disk limit between the Tevatron and 
LHC energies, see Fig. |2j Thus we expect that at the LHC the cross section will be less than 
that given by the DL parametrization, while the slope of elastic cross section, B e \, should be 
larger. The partons occurring in the proton wave function will be pushed away from the centre 
to the periphery. This suppression of the amplitude at low values of b automatically increases 
the elastic slope B e \, since B e \ oc R 2 = (b 2 ) where R is the interaction radius. 

For this reason the value of a' = 0.25 GeV~ 2 in (j5J) should not be considered as the slope 
of the bare Pomeron trajectory. Part of this slope is generated by the stronger absorptive 
corrections at smaller b. Indeed, the shrinkage of the diffractive cone observed in deep-inelastic 
scattering at HERA (say, in p-meson or in J/ip diffractive production) is smaller than that 
observed in pp-scattering. 

Actually to satisfy elastic unitarity it is sufficient to consider an eikonal model which sums 
up the exchanges of any number of Pomerons between the incoming protons. The model gives 
a result analogous to (J7J) 

A(b) =i(l -exp(-0(6)/2)) (9) 

where, up to a normalization factor, the opacity Q(b) plays the role of the scattering phase 5i. In 
the eikonal model, ([9]), it is the opacity Q oc s Qp_1 which is now described (or parametrised) by 
Pomeron exchange, rather than the final amplitude A. The situation is sketched symbolically 
in Fig. E](a). However this is not the whole story. The eikonal model does not include events 
with Large Rapidity Gaps (LRG), arising from the dissociation of one or both of the protons 
into high-mass systems, which are observed experimentally. We speak of single or double high- 
mass diffractive dissociation. We introduce high-mass dissociation in the next Section and give 
detailed discussion of these LRG events in Section 

Before we do this, let us recall that the interaction radius (i.e. the elastic slope) expected 
at the LHC will already exceed 1 fm - the distance where confinement may enter the game. 
Thus it is possible that confinement will stop the growth of the radius; that is, the growth of 
B e \. Another interesting possibility, proposed by V.V.Anisovich [10], is that the partons (the 
gluons and quarks) in the proton will start to form a few colourless clusters (like the nucleons in 
nucleus), rather than the normal homogeneous distribution. This may lead to a rich diffractive 

5 In general, an amplitude with 1 < ImA < 2 is not forbidden. The corresponding model, so-called U- 
matrix unitarisation, was discussed in [8]. However, such a model which asymptotically leads to a pure elastic 
interaction, without any particle production, does not look probable. 
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Figure 2: The impact parameter profiles of the imaginary part of the elastic amplitude for 
both the DL effective Pomeron pole [7] and the KMR multi-Pomeron [U] parametrizations, for 
collider energies y/s = 62.5, 1800, 14000 GeV relevant to the CERN-ISR, the Tevatron and 
the LHC respectively. 

dip structure in the differential elastic cross section da e \/dt. It is an additional argument to 
study experimentally the t behaviour of the elastic cross section at the LHC@. 

3 Diffractive dissociation 

So much for elastic scattering, which we may call elastic diffraction. Now we turn to inelastic 
diffraction, which is a consequence of the internal structure of the protons. Besides the pure 
elastic two-particle intermediate states shown in Fig. [3](a), there is the possibility of proton 
excitation, p — > N*, shown in the small sketch in Fig. [3]^b). At high energies, when the lifetime 
of the fluctuations of the fast proton is large, r ~ E/m 2 , the corresponding Fock states can 
be considered to be 'frozen'. Each constituent of the proton can undergo scattering and thus 
destroy the coherence of the fluctuations. As a consequence, the outgoing superposition of 
states will be different from the incident particle, so we will have inelastic, as well as elastic, 
diffraction. 

6 Reviews of the predictions for the total cross section, and the elastic scattering cross section over an extended 
\t\ interval, can be found in [11) . 
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Figure 3: (a) The single-channel eikonal description of elastic scattering; (b) the multichannel 
eikonal formula which allows for low-mass proton dissociations in terms of diffractive eigenstates 
\4>i), \4>k)', and (c) the inclusion of the multi-Pomeron-Pomeron diagrams which allow for high- 
mass dissociation. 



To discuss inelastic diffraction, it is convenient to follow Good and Walker [T2] , and to 
introduce states (pk which diagonalize the A matrix. These, so-called diffractive, eigenstates only 
undergo 'elastic' scattering. To account for the internal structure of the proton we, therefore, 
have to enlarge the set of intermediate states, from just the single elastic channel, and to 
introduce a multi-channel eikonal. The situation is pictured in Fig. E](b). 

What about proton dissociation into high-mass systems? At first sight, it appears that we 
may account for it by simply enlarging the number of diffractive eigenstates <pk- Even if this 
were practical, we would still face the problem of double counting when partons originating 
from the dissociation of the beam and the 'target' protons overlap in rapidities. Instead, high- 
mass dissociation is described in terms of so-called "enhanced" multi-Pomeron diagrams. The 
first, and simplest, is the triple- Pomeron diagram, shown in Fig. [I] and again in Fig. [3](c). In 
fact, high-mass dissociation is much larger than low-mass dissociation at very high energy, and 
it will play a central role in our discussion of soft interactions. For simplicity, therefore, we 
do not discuss low-mass dissociation further, but consider that it can be easily allowed for in 
terms of a multi-channel eikonal. So, from now on we should regard Q as an effective opacity 
embodying low-mass dissociation. 

Why does high-mass dissociation become so important at high collider energies? A simplified 
way to see this is to note that the cross section for dissociation of a proton into a high-mass 
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(M) system has the approximate fornfl 



f M 2 da SD dM 2 

with A = gzp/gN, where g 3 p is the triple- Pomeron coupling and is the coupling of the 
Pomeron to the proton, see Fig. [Q The Ins 'rapidity' factor comes from the integration 
J dM 2 /M 2 . Note that here we have used a e i ~ g% and M 2 da^jdM 2 ~ g%g%p ~ ^<7jV5 see 
Fig[H Eq. fflUj) is grossly oversimplified. For one thing, the dissociation cross section is subject 
to even greater absorptive (that is, multi-Pomeron) corrections than the elastic cross section. 
When these are taken into account it is found from analyses of the available high-energy data 
that 

X = g 3P /g N ~ 0.25, (11) 

see, for example, [9]. The important point is that the relative size of the a e \ and o"sd cross 
sections is governed, not simply by the parameter A, but rather by Alns ~ 1. For each fixed 
rapidity interval the probability of high-mass dissociation is relatively small, but the cumulative 
effect in the complete interaction amplitude is enhanced by the large phase available in rapidity 
at high collider energies. 



4 The structure of the Pomeron 

In the Regge approach, we have seen that high-energy soft interactions are driven by Pomeron 
exchange, together with its absorptive (multi-Pomeron) corrections. We may call this the "soft" 
Pomeron. Note that the s-channel asymptotic behaviour arising from t-channel Pomeron ex- 
change corresponds to an interaction radius which grows with energy, R 2 ~ Ins. The Pomeron 
must also describe multiparticle interactions. The sum of ladder diagrams of the type of 
Fig. Ufa) is the simplest multiparticle structure which reproduces the power-like s a behaviour of 
the Pomeron pole. Before the advent of QCD, ladder diagrams composed of t-channel mesons 
(dominantly spinless pion exchanged) were summed, and it proved difficult to generate an in- 
tercept a(0) > 1. 

Shortly after the discovery of QCD it was noticed [T3] that (colourless) two-gluon exchange 
has the properties of Pomeron exchange: vacuum quantum numbers, even signature and a 
singularity at j = 1. Due to the polarisation vectors of the spin-one particle (the gluon), the 
numerator of the propagator contains the g^ v tensor such that the 4-momenta of the incoming 
fast protons occur in the form p^g^p^ = (p a 'Pb) — s/2. Therefore the gluon-exchange amplitude 
contains an extra power of s as compared with spinless particle exchange. 

7 Here, for simplicity, we assume an essentially flat energy dependence, a ~ s e with elns < 1. The final 
equality in (fTO]) can be deduced by taking the ratio of the couplings indicated in the Regge expressions in Fig. [1] 

8 Pion-cxchange was considered to be the dominant mechanism due to the proximity of the pion pole, at 
t = mj, to the s-channel physical domain ,t < 0, due to the smallness of the pion mass. 
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Later, by applying the Regge approach to gluons, rather than to hadrons, it was possible, 
if the gluon transverse momenta k t were sufficiently large, to describe high energy (low x) 
interactions in terms of perturbative QCD. Then, on summing the leading contributions, where 
the smallness of the QCD coupling a s is compensated by large values of ml/a;, we build up 
the "hard" Pomeron. Below, we briefly trace this history and then apply these ideas to gain 
insight into high-energy soft interactions. 



4.1 Ladder structure of the Pomeron pole before QCD 

In terms of Regge theory, the high-energy amplitude is described by the exchange of a Pomeron 
(or a few Pomerons). As mentioned above, from a quantum field theory viewpoint, the Pomeron 
may be regarded as the sum of ladder-type diagrams, Fig. BJa). In other words, it corresponds 
to a sum of completely inelastic 2 — > n processes; that is, to the last term G; ne i = 1 — exp(— Q) in 
the unitarity equation ([6]). Now let us 'cut' the Pomeron; see Section [5721 for further discussion. 
That is, we treat the graph of Fig. BJa), not as a diagram for the amplitude, but rather as the 
diagram for the 2 — > n cross section 



a 



(2 -f n) = A*(2 -> n) ■ A(2 -> n), (12) 



see Fig. BJb). Then, we get the inelastic production of n particles homogeneously distributed 
in the available rapidity, y, interval. 

In the original 'soft' Regge approach it was assumed that the transverse momenta k t of 
these secondaries were limiteqj- The original idea, proposed by Amati, Fubini and Stanghellini 
in 1962 |15j . was to explain a large high-energy cross section, which did not decrease with 
energy, as a sequence of interactions in which the energy of a pair of adjacent particles (pions) is 
relatively small and lies in the resonance region where the interaction is strong. The integration 
corresponding to each cell of such a ladder is over the rapidity y and the transverse momentum 
kt, see Fig. H(a). If we denote the momentum transferred through the whole ladder by q t , then 
the contribution from a particular cell is a(q^)dy where 

( 2 ) Jl- [ fl-Vl 

a{qt> ~ 16tt 2 J (k 2 t +m?){{k t -q t y + m?y 1 ' 

and g is the coupling. The integration over the rapidities gives 

= T,^.flf ad vi = E^r = e " Y = (14 > 

n i n 

where the total rapidity spanned by the ladder Y = In s, and n\ reflects the identity of the cells 
or, equivalently, the ordering of the rapidities yi. 



9 Nowadays, with QCD, as mentioned above and discussed below, the Pomeron is described by a ladder built 
of the gluons. As a result, due to the dimensionless QCD coupling, there is no parameter to limit the value of 
kt, and BFKL equation predicts the 'diffusion' of the gluons in h\k t space [14] . 
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(a) (b) (c) (d) (e) 

Figure 4: Feynman diagrams corresponding to (a) one-Pomeron exchange, (b) the 2 — > n cross 
section, (c) two-Pomeron exchange, and (d,e) the triple-Pomeron contribution to high-mass M 
diffractive dissociation. The particles that are 'cut' by the vertical dashed lines in (b) and (c) 
are on-mass- shell. 



A convenient alternative way to obtain this result is to consider the evolution equation in 
rapidity 

^ = a A(y). (15) 
dy 

In such a form a plays a role similar to a DGLAP splitting function. 

Note that the value of the 'effective spin' a(qf ) depends on the momentum transfer through 
the ladder, —t = q 2 . The slope of the Pomeron trajectory a' = da{t)/dt at t = is controlled 
by the mean value of the transverse momenta a' oc g 2 / (k 2 ) or, in the case of ( |T3|) . by the mass 
m of the t-channel particle. 

However, with this procedure it turns out to be practically impossible to obtain a large 
intercept a(0) > 1 whilst keeping the t-channel exchanges to have masses m < 1 — 1.5 GeV. 
The corresponding couplings g are too small. On the other hand, if we include heavier states, 
then we go to larger kt] that is to small distances where it looks reasonable to work in terms of 
perturbative QCD. 



4.2 Ladder structure of the Pomeron pole after QCD 

With QCD, the Pomeron, at least at small distances, is described by the ladder of Fig. |4](a) 
built of gluons, rather than of spinless particles; with the coupling g 2 in (|T3|) replaced by the 
QCD coupling a s . The integral over k t now has a more complicated structure. It is convenient 
to consider the recursion relation between the amplitudes f n (x, k t ) with n and n — 1 cells 

fn(x,k t ) = ^ t% [^Kfakftfn-lVX), (16) 
71 J x X' J 71 

where we keep just the leading lnl/x' term in the x' integration. Note that now we have to 
account for the dependence of the amplitude on the transverse momenta, k t , of the t-channel 
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gluons. The recursion relation may be rewritten in differential form, as the BFKL evolution 
equation [T5] 

df N c a s 



d\nl/x n 



-K®f (17) 



so that / oc x A , where 



A = N c a s (K)/n = a(t) - 1, (18) 

where (K) is the leading eigenvalue calculated using equation ( !20l) below. Since we now have 
spin-one gluons in the ladder, the intercept a(0) = 1 already in the limit a s — > 0. 

In fact, if just the leading lnl/x contributions are summed then, explicitly, A = A = 
a s 41n2, where a s = 3a s /n [TB]. However, the next-to-leading logarithm (NLL) contributions 
are numerically rather large. It was found [TTj that, accounting these a s corrections, A ~ 
Ao(l — 6a s ). For the relevant, not too small, values of a s , a resummation is necessary. Since 
the NLL contribution is mainly of kinematic origin, it is possible, and looks reasonable, to 
resum all the major higher-order corrections. This leads to A ~ 0.3 [TB] in a wide region of k t . 

For simplicity, we consider just the forward amplitude^! with q t = 0. Then the elastic 
forward amplitude is 

A = i * I n& ( 19 ) 

J 7Tk t 

and, at LO, the BFKL kernel K acts as 

*<*,, K) 9 /(,, K ) = {/(,, K ) - k ^% ] ■ (20) 

Since in QCD we deal with a dimensionless coupling a s and a massless gluon, the mean value 
of kt in each cell is determined by the value of transverse momentum k' t in the previous cell. 
This is an important property of the BFKL equation - that is, the so-called diffusion in lnfc^. 
The value of L = Inkf in the current cell may differ from that in the previous cell by some 
quantity SL ~ 1. In other words, at each step of evolution, not only the impact parameter b 
can be changed by 5b ~ l/&f , but also the value of In k\ can be changed by 5L ~ 1. 



4.3 The transition from the 'hard' to the 'soft' Pomeron 

There are phenomenological hints that at large distances the "soft" Pomeron should have 
qualitatively similar structure as the "hard" (QCD) Pomeron. Indeed, first, no irregularity is 
observed in the HERA data in the transition region, Q 2 ~ 0.3-2 GeV 2 , between the 'soft' and 
'hard' interaction domains; the data are smooth throughout this region. Second, a small slope 
a' P ^0.1 GeV~ 2 of the Pomeron trajectory, is obtained in global analyses of all available soft 
high-energy data, after accounting for absorptive corrections and secondary Reggeon contribu- 
tions. This indicates that the typical values of k t inside the Pomeron amplitude are relatively 

10 For non-zero t the kernel K is a bit more complicated, but qualitatively similar in form. 
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large (a' oc 1/kf). Finally, recent 'soft' model data analyses [9j [19j [20] which account for the 
enhanced absorptive effects find an intercept of the initial, bare Pomeron A = ap(0) — 1 ~ 0.3 
close to the intercept of the BFKL Pomeron after the NLL corrections are resummed [IB] . Thus 
it looks reasonable to assume that in the soft domain we deal with the same perturbative QCD 
Pomeron; at least, there is a smooth transition from the soft to the hard Pomeron. 

In summary, the bare hard Pomeron, with a trajectory with intercept A = ap(0) — 1 ~ 0.3 
and small slope a', is subject to increasing absorptive effects as we go to smaller k t which allow 
it to smoothly match on to the attributes of the soft Pomeron. In the limited energy interval 
up to the Tevatron energy, some of these attributes (specifically those related to the elastic 
amplitude) can be mimicked or approximated by an effective Pomeron pole with trajectory 
intercept A = a P (0) - 1 ~ 0.08 and slope a' = 0.25 GeV~ 2 . 

5 Events with Large Rapidity Gaps 

For further discussion, we write the BFKL evolution equation (fTTj) in the oversimplified form 

*M = A A(y). (21) 

where A = ap(0) — 1 acts as the splitting function. This evolution equation is the analogue of 
( fl5!) . but, now, the possibility of multiparticle production leads to an additional power growth 
of the amplitude, A(t = 0) oc (s A )s, in comparison with the Born amplitude, Ab oc s. In other 
words, the evolution (!2T1) describes the development of the parton cascade from the 'beam' 
proton to the 'target' proton, and the s A growth of the inelastic cross section reflects the power 
growth of the parton density in the cascade. 

After one of the partons interacts with the target the coherency of the beam proton wave 
function is destroyed and a number of secondaries is produced. Note that in the ladder of 
Fig. IH(b) only one branch of the whole cascade is shown, which finally leads to a parton being 
absorbed by the target. The number N of secondaries produced is equal to the number of 
steps (partons) in this particular branch, shown in Fig. |4^b). According to fT2"T]) . where A plays 
the role of the particle density per unit of rapidity, we have N ~ A • Y. On the other hand, 
each of the partons continue to develop its own cascade. One example is shown in Fig. 11(b), 
where parton c develops its own branch. However, the branches which were not affected by the 
target conserve their coherence, and in the final state appear as a single parton. Thus we get 
the power growth of the amplitude, but only a logarithmic growth of the multiplicity N. The 
studies of coherence phenomena in a partonic cascade were pioneered by V.N. Gribov in 1972, 
see [21]. 

5.1 Large Rapidity Gap contribution to the total cross section 

The fact that the coherence of the wave function of the beam proton was destroyed after the 
interaction with the target, leads not only to inelastic high-multiplicity production, but via 
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unitarity relation (Q, also to elastic scattering. Elastic scattering is due to the absorption of an 
initial (coherent) component in some domain of the impact parameter b plane, and originates 
from the remaining part of the initial wave function which conserves the coherence between 
the partons. That is, elastic scattering occurs as the 'shadow' of the inelastic interaction. If 
the inelastic contribution Gi ne i is represented by the single ladder of Fig. BJb), then the elastic 
contribution may be drawn as a 'cut' between the two Pomerons of Fig. IH(c), see also Fig. [HI 

So much for elastic scattering of the incoming proton. However, an intermediate parton 
in a ladder may be scattered elastically as well. A unitarity equation, analogous to (El), may, 
and should be, written for any intermediate parton. The elastic scattering of an intermediate 
parton c can be represented by the diagram of Fig. H](d). Here, both branches of the cascade 
which start to evolve from parton c save their coherence. Therefore in the final state we will 
have a Large Rapidity Gap between the parton c and the target. Such a process is called the 
diffractive dissociation of the beam proton into a high-mass state M. It is described by the 
triple-Pomeron diagram of Fig. Hte), in which one Pomeron is split into two Pomerons, see also 
Fig. [TJ The probability r(= g^p) of this splitting within a unit rapidity interval is relatively 
small. First, due to the small parton c density A, where A arises from fl2T]) . Next, due to the 
parton-(target) proton cross section being smaller than the proton-proton cross section, since 
the incoming proton contains many partons. On the other hand, each intermediate parton 
in the ladder may generate a Pomeron splitting. Thus, the whole effect accumulated during 
the evolution (l2Tj) is enhanced by the parton multiplicity iV - in other words, by the size of 
available rapidity space Y, see also the discussion around equation (fit)]) . Therefore a triple- 
Pomeron diagram (like that in Fig. 11(e)) or diagrams with more complicated multi- Pomeron 
vertices (like those in Fig. [5(b) ) are called enhanced diagrams. 

Already at the Tevatron, the cross section of the events with a LRG, that is of single- 
and double-diffractive dissociation of one or both protons (integrated over the mass M of the 
final 'diffracted states'), is comparable to the elastic pp cross section. If we were to take the 
contribution of the diagram of Fig. H](e) literally and integrate over the rapidity of vertex r, 
then, at sufficiently high energy, the cross section of events with a LRG would exceed the whole 
inelastic G in ei contribution described by a single Pomeron, that is the contribution of Fig. H](b) 
(recall a P (0) > 1). 

The situation is even worse when we consider events containing a large number of LRG, 
arising, for example, from the multi- Pomeron diagram of Fig. E(a). For simplicity, let us first 
account for gaps of limited size, say 5y.i < y$. We denote the probability to form such a gap 
as 7. Then, integrating over the central position yi of each gap, we obtain for Y ^ y the 
contribution 



which grows as a power of s. If we neglect the limit 5yi < yo, and allow for the integration 
over the gap size Syi, then the cross section increases faster than any power of s, and clearly 
violates unitarity. This problem was discussed long ago in [22]. 
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Figure 5: A diagram for a process with several LRG centred at rapidities yf. (a) neglecting the 
gap survival factor S 2 , and (b) including the S 2 suppression caused by both eikonal (depicted 
here by the Pomeron line connecting the two protons) and enhanced rescattering involving 
intermediate partons. Plot (c) shows the expected long-range correlations R2{y a ,yb) of (f25j) . 
corresponding to diagram (a), with y a = y 1 , supplemented by eikonal screening. 

Two main scenarios were proposed to prevent such a growth and to restore unitarity. First, 
the so-called weak coupling scenario assumes that the triple-Pomeron vertex r oc t (and corre- 
spondingly the value of 7) vanishes as t — > [23J. Since the mean value of momentum transfer 
it) decreases as 1/ In s (the shrinkage of the diffractive cone) the decrease of r compensates 
the increase of the available rapidity interval. However this hypothesis is not supported by 
experimenj^l. 

In the favoured strong coupling scenario |25j, the multi-gap cross section is suppressed by 
a small 'gap survival probability' S 2 which decreases with energy. In this case the Pomeron 
coupling r — ► const as t — > 0, while the small probability of LRG production is due to the large 
probability that the gaps are populated by secondaries produced in additional soft rescattering 
interactions between the protons (and also the intermediate partons). 

In the simple eikonal model, 0, the probability, not to have an additional inelastic interac- 
tion which will populate the gap, is given by a factor S 2 = exp(— Q). Recall that in the eikonal 
model Ginei = 1 — exp(— Q). However, this factor accounts for the rescattering of the incoming 
protons only. Since, now, we have also rescattering between the intermediate partons we have 
to consider more complicated multi- Pomeron diagrams, like those shown in Fig.3(c). 

11 Still the weak coupling possibility is not rejected completely and it will be important to study high-mass 
dissociation at the LHC in the low t region to confirm this conclusion; see [24] . 
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Figure 6: Different processes described by cutting the two-Pomeron exchange diagram. The 
four diagrams correspond, from left to right, to the contribution to the imaginary part of the 
elastic amplitude; to the cross section o~ with zero multiplicity in the central region; to the 
process with single Pomeron multiplicity (01); and, finally, with double Pomeron multiplicity 



5.2 Multi-Pomeron diagrams and the AGK cutting rules 

We emphasize that each multi-Pomeron exchange diagram describes simultaneously a few dif- 
ferent processes. We have already seen that the first ladder diagram gives, on the one hand, 
the elastic pp scattering amplitude, Fig. H^a), while, on the other hand, it may be considered 
as the 'cut' diagram for the cross section for multiparticle production, Fig. IU(b). In the latter 
case cutting the Pomeron of Fig. S](a), that is taking the discontinuity, discv4=2Iim4, gives the 
contribution to multiparticle production, Gi ne i = 21mA 

Now consider the two-Pomeron exchange diagram of Fig. IU(c). Strictly speaking, the dia- 
gram should be drawn more precisely, as has been done in Fig. We have the possibility to 
cut both Pomerons simultaneously, 0(d), to cut only one Pomeron, 0(c), or to cut between the 
two Pomerons, 0(b). The sum of all cuts gives the total contribution of two-Pomeron exchange 
to the imaginary part of the elastic amplitude. This contribution is negative and describes 
the absorptive correction to the one-Pomeron amplitude. On the other hand, the elastic cross 
section corresponding to Fig. 0(b) is positive; as is the cross section of the events with a 'dou- 
bled' particle density, which we obtain by cutting both Pomerons as in Fig. 0(d). The only 
negative contribution is the absorptive correction to the events with a single particle density 
(i.e. the correction to single- Pomeron exchange). The famous AGK cutting rules gives the 
relation between the different subprocesses originating from the same Reggeon diagram. For 
the respective cut diagrams of Fig. 0, the relation is 

(7 : o x : 0-2 = 1 : -4:2, (22) 

while the whole contribution to the elastic amplitude isl — 4 + 2 = — 1. 

Analogous relations hold for more complicated diagrams with many Pomerons. Note that, 
according to the AGK rules, the multi-Pomeron diagrams do not change the single particle 
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inclusive cross section da/dy. Indeed, the particle density arising from diagramed) is twice 
larger than that from^c), and these two contributions cancel each other. The same is true for 
more complicated multi-Pomeron diagrams. 

If we start with the eikonal model, then the sum of all the multi-Pomeron exchange diagrams 
gives the total probability of the interaction at fixed b, 

atot = 2(1 - e-^) = 2f>ir 1 ^^ = f>(&), (23) 

n=l n=0 

where Q(b)/2 is given by one-Pomeron exchange and 

on 

a = a el = (1 - e^/ 2 ) 2 , a n = ^e~ n for n > 1. (24) 

n! 

Here, <7 (b) is the probability of the elastic interaction; and a n (b) is the probability to produce 
a particle density n times larger than that arising from one-Pomeron exchange, that is the 
probability to cut n Pomerons in the whole amplitude. It is easy to see that 

oo oo ^ n 

a inel = °n = ~[ e ~ n = l-e- n = G inel (25) 

n=l n=l 

is the whole inelastic contribution at a fixed b. The penultimate expression for cr inel (6) rep- 
resents the whole probability, 1, minus the probability e~ n to have no inelastic interaction, 
whereas in the previous expression for a ine \ each term Q n /n\ represents the probability of n 
inelastic interactions (where n\ accounts for the identity of the interactions) multiplied by e~ n 
which guarantees that there are no further inelastic interactions. The multi-Pomeron exchange 
structure of processes of ( I2"3l -( l2~5l) is sketched in Fig. [7J 

Similar rules, with the same combinatorics, are valid in the presence of enhanced screening. 
At each multi-Pomeron vertex, we can cut one or more Pomerons, or cut between the Pomerons 
placing the uncut Pomerons either to the left or to the right of the cut. For illustration, we 
discuss the contributions arising from cutting a triple-Pomeron diagram which is screened by 
a single eikonal Pomeron, as shown in Fig. [8](a). The Pomeron ladders corresponding to this 
diagram, Fig. \EHjo), can be cut in 8 different ways. First, we have the elastic cxo cut between 
all the ladders. Next, the cut of the eikonal Pomeron ladder gives the usual single- Pomeron 
multiplicity. Then, we have three possibilities to cut the triple-Pomeron graph; one with single 
multiplicity in the whole rapidity interval, one with a gap in the upper part, one with double 
multiplicity in the upper part. Finally, we have the same three, but now simultaneously cutting 
the eikonal Pomeron. Thus, the different cuts, which lead to 8 different process, give the particle 
multiplicities as a function of rapidity that are sketched in Fig. E^c). The AGK cutting rules 
may be used to determine the relative probabilities of these processes; the result is given below 
the multiplicity plots. For the example of Fig. [8]we have to use relation (1221) twice. For instance, 
the probability —8 comes from 2 x —4, whereas —4 comes from —2 x 2 (—2 corresponds to 
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Otot = XI 



n=l 



0/2 



n 



ni,n2=l 



ni,ri2=0 




Figure 7: The eikonal structure of the total interaction: a tot (b) = cr e i(&) + cr inel (6) = a (b) + 
X] n cr n (6). The dashed vertical lines represent the cut leading to the elastic contribution, a , 
and the cuts of the n Pomerons in cr„. The circle is to indicate that all n Pomerons are cut. 



the cut of the triple-Pomeron diagram without cutting the single Pomeron; note that in (1221) 
—4 = —2 — 2, where the two terms correspond to cutting one or the other Pomeron in Figj6]). 

Moreover, AGK rules show that, in general, these factors do not depend on exactly how 
the multi-Pomeron vertices are cut. To be specific, the probability of a particular subprocess 
is given just by the combinatorial factors to choose n cut Pomerons from m = n\ + n + n<i 
Pomerons of the original amplitude. A factor 2 comes from discA=2ImA of each cut Pomeron, 
and another factor 2 comes from the possibility to place any uncut Pomeron to the right or left 
of the cut. 



5.3 Long-range rapidity correlations 

Note that the eikonal model predicts a long-range correlation between the secondaries 
produced in different rapidity intervals. Indeed, we have possibility to cut any number of 
Pomerons. Cutting n Pomerons we get an event with multiplicity n times larger than that 
generated by one Pomeron. The probability to observe a particle from a diagram where n 
Pomerons are cut is n times larger than that from the diagram with only one cut Pomeron. The 
observation of a particle at rapidity y a , say, has the effect of enlarging the relative contribution 
of diagrams with a larger number of cut Pomerons. For this reason the probability to observe 
another particle at quite a different rapidity y\, becomes larger as well. This can be observed 



17 




Figure 8: (a) A triple- Pomeron diagram screened by one eikonal Pomeron. (b) The gluon ladder 
structure of the diagram. The dashed line and circle imply that there are several different ways 
(in fact, 8) to cut this multi- Pomeron diagram, (c) The particle density dN/dy, as a function of 
rapidity y, expected for the processes obtained from these 8 different cuts, together with their 
relative probabilities, which may be obtained from the AGK cutting rules. 



experimentally via the ratio of inclusive cross sections 

(T inc id 2 a/dy a dy b _ d 2 N/ dy a dy b _ 

2 1 : (da/dy a )(da/dy b ) ~ (dN/dy a )(dN/dy b ) ' 1 j 

where dN/dy = (l/a ine {)da/dy is the particle density. 

Without multi- Pomeron effects the value of R 2 exceeds zero only when the two particles are 
close to each other, that is when the separation \y a — yb\ ~ 1 is not large. Such short-range 
correlations arise from resonance or jet production. However, multi-Pomeron exchange leads 
to a long-range correlation, R 2 > 0, even for a large rapidity difference between the particles, 
\Va ~ Ub\ ~ Y. In fact, such a long-range correlation, with R 2 ~ 0.2 for \y a — y b \ ~ 5, has even 
already been observed in the old CERN-ISR data [27]. 

In the eikonal case, the value of R2 may be calculated explicitly from the Poisson distribution 
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of a n given in (1241) . At fixed b it gives 



R-2 



- 1. 



(27) 



For a low optical density (Q <C 1) the correlation arises mainly from the two-Pomeron exchange 
diagram, and increases as R2 = Q/2. At very large Q, the Poisson distribution becomes narrow 
in the number n of cut Pomerons, and R2 decreases as 1/Q. The maximum value is R2 — 0.3 



It would be informative to study the Pomeron loop structure by observing long-range cor- 
relations at the LHC energy. In the pure eikonal case, the correlation should be the same 
throughout the whole available rapidity interval (except, of course, near the edges of rapidity 
interval, where multiparticle production is affected by energy conservation constraints). When 
we have enhanced diagrams, that is Pomeron loops which produce LRG which occupy only 
a part Sy of the rapidity space, the corresponding correlation takes place only within this 5y 
interval. Thus the existence of Pomeron loops may be revealed both as LRG events (from the 
cut with zero multiplicity inside the Sy interval), and as long-range R2 correlations in the 8yi 
intervals (from the cuts generating large multiplicities of secondaries). The 'range' of the cor- 
relation reflects the size of the Pomeron loop in rapidity space. An example of this behaviour 
of R2 is shown in Fig. [5]^c) 

5.4 The decreasing cross section 

Note that the simplest scenario to restore the unitarity at very high energies is not yet com- 
pletely rejected. We mean the possibility that asymptotically the Pomeron intercept becomes 
less than 1, «p(0) < 1, and at very high energies the total cross section starts to decrease 
with energy. Such a behaviour is expected in a theory with only the triple-Pomeron coupling, 
and which neglects the more complicated multi-Pomeron vertices, such as the 2^2 Pomeron 
coupling and so on |28j. 

Recall that the two-Pomeron loop gives a positive probability for LRG production but its 
contribution to the amplitude is negative, see Fig. [6j If we now sum up the diagrams with 
many two-Pomeron loop insertions (Fig. [5](a)) then we 'renormalise' the Pomeron propagator. 
In this way we obtain the 'dressed' propagator in terms of a Schwinger- Dyson equation. The 
intercept of the final 'dressed' propagator is Op(0) = 1 — e=l + A — 'loop'. With increasing 
energy, the available rapidity interval increases, with the result that the loop renormalisation 
grows and hence the effective intercept 1 — e decreases. Asymptotically, the Schwinger-Dyson 
equation gives 



for n 



1.8. 



e = 




(28) 
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where r is the triple- Pomeron coupling. The loop contribution, 'loop' = r 2 /2e, is calculated 
using the dressed Pomeron propagator, and the factor 1/e comes from the integral over the 
loop size 5y^ 

For the reasonable values r ~ 0.2 and A ~ 0.3, the naive estimate based on (T2"8"j) gives 
e ~ 0.06. That is ap Sympt (0) ~ 0.94. This means that up to a rather large energy (Y = Ins < 
1/e ~ 18) the total cross section will grow, but after this, that is just in the LHC region, it will 
start to decrease. For multiparticle production, such a regime will reveal itself as long-range 
rapidity fluctuations of the multiplicity arising from two-Pomeron loops which may produce 
either LRG or a double multiplicity of secondaries in the rapidity domain occupied by the loop. 



6 Saturation of the kt distribution 

Let us return to the simplified eikonal model (Q, written for the partons in the incoming 
protons. We may ask the question - "what happens to an incoming fast quark?" On the one 
hand, the factor exp(— Q) describes its absorption. After the inelastic collision the incoming 
quark disappears from the initial wave function. On the other hand, the quark cannot just 
disappear 3 ]. It has baryon charge, etc. Disappearance from the incoming beam (wave function) 
actually means migration - after the collision the momentum of quark is changed. In the 
leading logarithmic approximation, the quark mainly changes its transverse momentum. Each 
new Pomeron-quark coupling may be considered as a new elastic scattering and it is known 
that after n scatterings the mean transverse momentum squared, (k^ ), increases n times. 

Thus, when we say that unitarity stops the growth of the parton density via an increase of 
the absorptive correction, actually we mean that the extra partons produced by the cascade are 
pushed out of the previous domain in b, kt configuration space into a region of larger kt and b. 
Accounting for the enhanced diagrams we get the same phenomena, not only for the incoming 
fast partons (as in the eikonal model), but for any intermediate partons as well. Thus finally 
we expect that the particle density will reach saturation at low k t , but will continue to grow 
with energy at larger k t . Therefore, at higher energies a larger number of minijets with larger 
(k t ) will be produced. 

Due to saturation in the low k t region at very high energies, the majority of partons will have 
rather large k t . We already know that the diffusion in In kt takes place in the perturbative QCD 
domain [14J, and to provide a smooth matching to the larger k t domain we need to include some 
elements of a similar diffusion in the description of low (and intermediate) transverse momenta. 

12 Strictly speaking, there should also be a Schwinger-Dyson equation giving the dressed triple-Pomeron 
coupling. However, in this oversimplified estimate, we have truncated the system of Schwinger-Dyson equations 
at the first step. 

13 For the gluon the situation is more complicated. Some of gluons can disappear via fusion, gg — ► g. 
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7 Information from double- Pomeron-exchange reactions 



At the Tevatron and the LHC, there is a possibility to study just Pomeron-Pomeron collisions 
by selecting events with large rapidity gaps on either side of some centrally produced hadronic 
system of mass M. This can be done, either, by detecting the two forward outgoing protons 
which each carry a large fraction, xl : of their incoming momenta0, say xl > 0.96 , or, just 
by using the "gap-hadron-gap" trigger. Since the interaction across each gap is described by 
Pomeron exchange such events may be treated as Pomeron-Pomeron collisions. Usually it is 
called a Double-Pomeron-Exchange (DPE) process. 

In a DPE process we may hope to observe the properties of the Pomeron caused by the 
fact that (i) the Pomeron consists mainly of gluons and (ii) the transverse size of the Pomeron 
is small, that is the typical transverse momenta, kt, in the diagrams which describe Pomeron 
exchange are relatively large. Indeed, the value of k t increases due to BFKL diffusion and 
screening effects (as explained in the previous Section). Another indication in favour of large 
kt is the rather small value of the slope a' oc 1/ (k 2 ) of the Pomeron trajectory, see Section fl~3l 
and the discussion after eq. (|T5|) . 

For these reasons, DPE events have several characteristic features. First, since the Pomeron 
consists of mainly of gluons, we expect there to be a larger fraction of the glueballs and/or rj', 77 
mesons in the central system M 2 , than that in pp-interactions. Recall that SU(3) iso-singlet 
7/, 77 mesons contain a sizeable gluon component. Moreover, because DPE has no "incoming" 
valence quarks, we expect the baryon/hyperon content to be smaller than in pp interactions. 

Further information can be obtained by comparing DPE processes with pp interactions at 
the lower energy yJ~Sp~p~ = M. For example, the multiplicity of secondaries is expected to be 
close to that in the pp case0. Strictly speaking, we expect a bit higher multiplicity since in 
Pomeron-Pomeron collisions the initial energy goes mainly to mesons, while the fraction of 
baryons is smaller than in pp interactions. 

There are several consequences of the small size of the Pomeron. First, the transverse 
momenta of secondaries are expected to be a bit larger for DPE than that in the pp case at 
s pp = M 2 , but smaller (due to the kinematic constraints) than in the original pure inelastic 
process which takes up the whole initial pp energy */s. The high Et dijets produced by DPE 
will be dominantly gluon jets. Next, the relatively small size of the Pomeron should result in a 
smaller interaction radius, as measured via Bose-Einstein correlations, than that found in the 
usual inelastic pp collision. Finally, the probability of double parton scattering, in particular 
the probability to observe two pairs of high Et jets, will be larger than that in a pp collision 
at Spp = M 2 . 

14 The mass of the centrally produced hadronic system is given by M 2 = s(l — ££,i)(l — Xlz) and the size of 
the gaps in rapidity are Syt = — ln(l — xn). 

15 We cannot compare with e + e~ interactions, since there the multiplicity, and other characteristics of the 
final distribution, are driven by specific double loagarithms which are different from those in pp collisions |29j . 
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Accounting for the anticipated LHC luminosity, the rate of DPE events is rather large. For 
example, from Ref. [9], we see that the DPE cross section integrated over the 0.002 < £ < 0.02 
interval, corresponding to the central production of a system of mass M ~ 100 GeV, is about 
10 \ib. Alternatively, the expected DPE cross section is ii&da / d^d^ ~ 1 — 5 \xb [9j; here 
£i = 1 — xu- So it should be possible to study the properties of DPE events in detail. 

8 Conclusions 

We should discuss how the saturation of particle densities can be consistent with confinement. 
The growth of the transverse momenta k t tames the increase of the interaction radius with 
energy. One possibility is that the simultaneous action of confinement and the growth of k t will 
finally lead to the saturation of the interaction radius; that is to the saturation of the elastic 
slope, B e \ — * const as s — ► oo. 

However the idea, discussed at the end of Section [2J that when the interaction radius becomes 
large confinement leads to the formation of colourless clusters, allows an alternative possibility. 
In a large, but finite, interval of rapidity the evolution is driven by the gluon ladder. In this 
case the impact parameter is almost frozen; that is, the variation of the interaction radius 5b is 
small. The system then emits a pion (or another meson) and this 'colourless cluster' provides 
a noticeable increase of the impact parameter b. Following this, we again have gluons in some 
interval of evolution, and so on, as shown in Fig. [5Ja). The hard QCD interactions will be 
responsible for the growth of kt, while the part of the evolution described in terms of colourless 
mesons will be responsible for the increase in the interaction radius, that is of B e \, due to the 
diffusion in b plane see Fig. Mjo). In such a 'mixed' picture we expect lower k t in peripheral 
collisions (since the available rapidity interval was used to enlarge the value of b) and maximal 
kt in central collisions. 

It may be informative to study multiparticle production at the LHC for individual events, 
as well as in terms of inclusive cross sections. At the high LHC energy the multiplicity of 
secondaries is large. Thus it might be possible to observe the fluctuation of particle densities in 
individual events like those shown in Fig. [SJ or the effects arising from a mixture of intervals of 
evolution via colourless meson and coloured gluon exchanges, like those shown in Fig. [9j In the 
latter case, the rapidity correlation length (measured via R2 of (126]) ) should be rather small for 
peripheral collisions since the main part of the available rapidity interval was used for meson 
exchange (to enlarge b) and so only small intervals are left for gluon exchange where we expect 
a large probability to have several 'cut' Pomerons. 

It is appropriate to say a few final words about models of high-energy soft interactions. We 
emphasize that the aim of experiment is not to reject or to confirm one or another model by 
saying that the data prefer Monte Carlo version X or Y. Rather, the objective should be to 
isolate, and to study, the main qualitative features of the interaction. What, therefore, are the 
requirements of a realistic model of high-energy soft interactions? To describe all the qualitative 
features discussed in this paper, it is clear that a realistic model should 
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Figure 9: The effect of the presence of colourless clusters (meson exchange) on the impact 
parameter of the produced particles as a function of rapidity. 



• contain, not only the eikonal but also the enhanced multi-Pomeron contributions; 

• include, not just the triple-Pomeron vertex, but more complicated multi-Pomeron vertices 
(otherwise we will get an asymptotically decreasing cross section); 

• allow for 'diffusion' of the partons both in impact parameter, b, and in transverse mo- 
menta, or, more precisely, In k t space. 

A recent attempt at building such a multi-Pomeron model, tuned to describe all the available 
data for high-energy soft interactions, is presented in [9]. For convenience, we call this the 
KMR model. 

Finally, can the characteristic features of the multi-Pomeron description of soft interactions 
be observed at the high Tevatron and LHC energies? We list some below. 

• The multi-Pomeron absorptive effects tame the power growth of the pp total cross section 
leading to values smaller than predicted before. For example, the KMR model prefers a 
value close to the lower of the two Tevatron measurements, and a value of about 90 mb 
at the LHC energy of ^ = 14 TeV, see Fig. ITDTa). At 100 TeV the prediction of 108 mb 
is at the lower limit of cosmic ray expectations. 

• Simultaneously, these absorptive effects shrink the elastic differential cross section peak 
as the energy increases through the LHC range more than the expectations arising from 
a naive effective Pomeron-pole of a Donnachie-Landshoff type of parametrization [TJ. 
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The cross section for proton dissociation into high-mass systems should grow with energy 
and be comparable to the elastic scattering cross section. The reason is that although, for 
each fixed rapidity interval, the probability of high-mass dissociation is relatively small, 
the overall effect is enhanced by the large phase space that is available in rapidity at high 
collider energies. The predicted values of <r e i and <tsd for the KMR multi-Pomeron model 
are shown in Figs. ITOTb.c) for energies in the Tevatron - LHC energy range. 

The main growth in multiplicity, as we go from Tevatron to LHC energies, is due to the 
small size ('QCD') Pomeron component, which produces particles with typically pt ~ 5 
GeV. There is essentially no growth in multiplicity at small p t . This simply confirms 
the trend that has been observed through the CERN-ISR to Tevatron energy range, see 
the data points in Fig. [TTJ In other words, starting with the same Pomeron intercept 
(A ~ 0.3 of Section |4~2]) the contribution of the large-size component after the absorptive 
correction becomes practically flat in energy (~ s 08 ), while the small-size contribution, 
which is much less affected by the absorption, continues to grow with energy (~ s ' 3 ). 
Such behaviour is consistent with the experiment (see Fig. [TTT) where the density of low 
p t secondaries is practically saturated while probability to produce a hadron with a large 
(say, more than 5 GeV) transverse momentum grows with the initial energy. 

Multi-Pomeron exchange diagrams have a characteristic pattern of 'enhanced' multiplic- 
ities of secondary particles, as well as large rapidity gaps, governed by the AGK cutting 
rules. These lead to long-range rapidity correlations, for example, which may be observed 
at the LHC via R 2 of (1251). 

The observed rate of Large Rapidity Gap (LRG) events in various processes at the LHC 
will be particularly informative. The survival probability of a LRG may be calculated in 
terms of a multi-Pomeron exchange model of soft interactions. It is necessary to calculate 
the survival to both eikonal soft rescattering (between the colliding protons), and enhanced 
rescattering (involving soft interactions with intermediate partons with different k t ). To 
calculate the latter we need to include the k t dependence of Pomeron exchange, see [9] - 
this was the third requirement listed for a realistic model. 

A topical example of a LRG process is the exclusive production of a heavy mass system 
A, that is pp — > p + A + p where the + signs denote LRGs. This process, with A =Higgs, 
is a novel and promising way to study the Higgs sector at the LHC, which gives a strong 
motivation for the addition of forward proton detectors to enhance the discovery and 
physics potential of the ATLAS and CMS detectors at the LHC [H [2] . Already, at the 
Tevatron, the CDF collaboration have observed exclusive processes with A = 77, dijet and 
Xc with rates consistent with the estimates of the LRG survival probabilities predicted 
by the multi-Pomeron model of soft interactions, see PQ and references therein. More 
data on exclusive processes, from both the Tevatron and, especially, the LHC, will be 
illuminating. 
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Figure 10: The energy dependence of the total (a), elastic and diffractive dissociation (b) 
pp cross sections, and the cross sections of dissociation to a fixed M 2 = £s state (c); (d) the 
parton multiplicity (solid lines) and the number of 'colour tubes' (dashed) produced by Pomeron 
components of different size. The figure is taken from the KMR model of [9]. 
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Figure 11: The plot is from Ref. [30J. The horizontal lines, which are superimposed, are the 
predictions of [9] at the Tevatron and LHC energies; the three p t ranges correspond to the 
large-, intermediate- and small-size components of the Pomeron, see also Fig. [LW d) 
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